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Abstract 

We analyze the onset of classical behaviour after a second-order phase tran- 
sition by considering a scalar field theory in which the system-field inter- 
acts with its environment, represented both by further fields and by its own 
short-wavelength modes. Within our approximations we see that the long- 
wavelength modes have become classical by the time that the transition has 
been first implemented (the spinodal time). 
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Cosmology and particle physics suggest strongly that phase transitions have occurred in 
the early Universe, in particular at the grand unified and electroweak scales 

An analysis of phase transitions in quantum field theory that takes the non-equilibrium 
nature of the dynamics into account from first principles is very difficult, and has only be- 
gun to be addressed. In particular, the naive picture of a classical order parameter (infiaton 
or Higgs) field rolling down an adiabatic effective potential, that was once a mainstay of 
cosmological field theory modelling, has been shown to be suspect . Alternatively, the sug- 
gestion by Kibble ^ that, while a non- adiabatic approach is crucial, causality alone can set 
saturated classical bounds on time and distance scales during a transition, has been shown to 
be only partly true. The issue of how a quantum system evolves into the classical theory has 
been addressed in Refs. 0,3. For some models, it has been shown that classicality emerges 
as a consequence of profuse particle creation, whereby a non-perturbatively large occupation 
number of long-wavelength particles produces, on average, a diagonal density matrix. This 
dephasing effect occurs at late times. Here we will consider a model of an explicitly open 
system, in which classicality is an early time event, induced by the environment. The result 
is completely different from that of Refs. [^,0 . 

There are several time scales which are relevant for the description of the onset of a 
transition. If the quench is fast, the initial stages of a scalar transition can be described 
by a free field theory with inverted potential, (mass)^ < 0. This description is valid until 
the field wave functional explores the ground states of the potential at the spinodal time 
tsp. Specifically, the field ordering after the transition is due to the growth in amplitude 
of its unstable long-wavelength modes. In consequence, the short-wavelength modes of the 
field, together with all the other fields Xa with which the inescapably interacts, form an 
environment whose coarse-graining makes the system-field classical ||^. As a result, there 
is an additional time scale associated with the environment, the decoherence time 
Once t > to the order parameter becomes a classical entity 0. 

In this letter we consider a simple model of a scalar order-parameter field 0, whose Z2 
symmetry is broken by a double-well potential. Specifically, we take the simplest classical 
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action with scalar environmental fields Xa (yU^i'^^a > 0) 



(1) 



^ f (1 1 ] 

5env[Xa] = J d^X j-^^Xa^^Xa - ^mlxlj , 

As we have already observed, the separation in (|I]) is not yet the separation between the 
ultimately classical 'system' and its 'environment' since the short- wavelength modes of the 
0-field never become classical, and must be treated as part of the decohering environment. 
However, to demonstrate how an environment renders the order-parameter field classical we 
consider the simpler case in which the environment is taken to be composed of the N fields 
Xa alone. Since the effect of environmental sources is additive, at our level of approximation, 
the Xa fields alone give us an upper bound on to- 

For weak couplings A, (?a ^ 1 and comparable masses ma ~ ;U we shall find that t_D is 
shorter than the spinodal time t^p, defined as the time for which ^ rf = the 
ground states of the system. In consequence, by the time that the field is ordered it can be 
taken to be classical. 

The model has a continuous transition at a temperature T^. The environmental fields Xa 
reduce Tc and, in order that T"^ = /i^/12(A + Z^fi'a) ^ /W^, we must take A + X^fl'a ^ 1- For 
order of magnitude estimations it is sufficient to take identical (7a = 9/y/N- We will also 
assume g — X. For one- loop consistency in our subsequent calculation we assume that 3> 
1. The effect of many weakly coupled environmental fields is that they act on the system field 
without it being able to act back on them. For example, the dominant hard loop contribution 
of the 0-field to the Xa thermal masses is 6m'^ = 0{gT^ /^/N) = 0{jj,'^/N) <^ jj,'^. Similarly, 
the two-loop (setting sun) diagram, which is the first to contribute to the discontinuity of 
the x-field propagator, is of magnitude g'^T^/N = 0{giJ? /N^^"^) <^ 5m^, in turn. That is, 
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the effect of the thermal bath on the propagation of the environmental x-fields is ignorable. 
We stress that this is not a Hartree or large-N approximation of the type that, to date, has 
been the major way to proceed for a closed system. 

We shall assume that the initial states of the system and environment are both thermal, 
at a temperature Tq > T^. We then imagine a change in the global environment (e.g. 
expansion in the early universe) that can be characterised by a change in temperature from 
To to Tf < Tc. The relevant object is the reduced density matrix Pi[(f>^, (f)~, t] = (0+|pr(t)|0~). 
It describes the evolution of the system under the influence of the environment, defined by 



where pf^"*", Xa > 0"; Xa > ^] — {^~^xt\pi't)\'P~Xa) the full density matrix. The environment 
will have the effect of making the system effectively classical once pT:{t) is, approximately, 
diagonal in the field configuration basis. Quantum interference can then be ignored and we 
obtain a classical probability distribution from the diagonal part of Pr{t), or equivalently, 
by means of the reduced Wigner functional, which is positive definite after the decoherence 
time (see Ref. for an explicit demonstration of this for a toy model). This behaviour is 
essentially different from the dephasing effects Our onset of classical behaviour is an 
early-time event which, beneficially, allows us to use perturbation theory. 

Assuming that the initial full density matrix can be factorised, the temporal evolution 
of the reduced one is given by 







where Jr is the reduced evolution operator 



J (h. J 6. 



(2) 



The Feynman- Vernon |jTT| influence functional F [0+ , ] is defined as 




X exp {z{S[xt] + 5int[0+,Xa+] " ^{^[Xa1 + 5intr,Xa1})- 
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Given our thermal initial conditions it is not the case that the full density matrix has (j) 
and X fields uncorrelated initially, since it is the interactions between them that leads to the 
restoration of symmetry at high temperatures. Rather, on incorporating the hard thermal 
loop 'tadpole' diagrams of the x (and </>) fields in the mass term leads to the effective 
action for (p quasiparticles, 

where m^(To) = — — Tq /T^) > 0. As a result, we can take an initial factorised density 
matrix of the form p[Tq] = p^[To]p^[Tq], where p</,[To] is determined by the quadratic part of 
SsysM and p^[To] by S[x^]. 

Provided the change in temperature is not too slow the exponential instabilities of the 
0-field grow so fast that the field has populated the degenerate vacua well before the temper- 
ature has dropped to Tf. As Tc has no particular significance for the environment field, for 
these early times we can keep the temperature of the environment fixed at = Tq = 0{Tc). 
Since it is the system-field field whose behaviour changes dramatically on taking through 
Tc, we adopt an instantaneous quench for from Tq to Tf = at time t = 0, in which m'^{T) 
changes sign and magnitude instantly, concluding with the value m'^it) = — yU^, t > 0. 

Beginning from this initial distribution, peaked around = 0, we follow the evolution 
of the system, with Hamiltonian determined from (|T|). From the infiuence functional we 
define the infiuence action 5y4[0+, 0^] by T[0^, 0^] = exp i6A[(f)'^, 0^]. After further defining 
A = |(0+^ — 0~^) and S = |(0"''^ + 0~^), the real and imaginary parts of the infiuence action 
are, in the one loop (two vertices) and large approximations^, 

Re6A =^ld''xl d% A{x)K^ix - y)S(y), 
ImSA =-^l d^'x I d^y Aix) N^ix , y) Aiy) , 



^The large-N approximation singles out the two-vertex loop in the one-loop perturbative 
approximation . 
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where -ft'q(x — y) = \mG\j^{x,y)0{y^ — x°) is the dissipation kernel and N^^i^x — y) = 
ReG^^(x,y) is the noise (diffusion) Icernel. G++ is the relevant closed-time-path correlator 
of the X" field at temperature Tq. 

The first step in the evaluation of the master equation is the calculation of the density 
matrix propagator from Eq.(0). In order to estimate the functional integration which 
defines the reduced propagator, we perform a saddle point approximation 

^r[0f^,0r,t|0+,0r,^o] ~ expiA[0^,0^], 

where 74[0+,0^] = S'[0+] — S[(j)~] + 5y4[0+, </)^], and 0^ is the solution of the equation of 
motion = with boundary conditions 0^(to) = 0f and = (t>f- It is very 

difficult to solve this equation analytically. For simplicity, we assume that the system-field 
contains only one Fourier mode with k = k^. We are motivated in this by the observation 
T0| , p!2[| that the exponentially growing long-wavelengths increasingly bunch about a wave- 



number /cq < which diminishes with time initially as /cq = 0{fi'^/t). 

The classical solution is of the form (j)ci{x,s) = f{s,t)cos{kQ.x) where /(0,t) = 0i and 
f{t,t) = 0f. Qualitatively, f{s,t) grows exponentially with s for t < tgp, and oscillates for 
tsp < s < t when t > t^p- We shall therefore approximate it, for t < tgp as 

^\ , sinh(u;os) , sinh[uJo{t - s)] 

f{s,t) = (1)1 . + 0i . , . — , 

smh[uJot) smh(ci;otj 

where uJq = fi"^ — k"^. In order to obtain the master equation we must compute the final time 
derivative of the propagator Jr, and after that eliminate the dependence on the initial field 
configurations 0f coming from the classical solutions 0^ (see Ref. [||). 

To determine the onset of classical behaviour it is sufficient to calculate just the correction 
to the normal unitary evolution coming from the noise kernel. For clarity we drop the suffix 
f on the final state fields. If A = {(f)'^'^ — 0~^)/2 for the final field configurations, then the 
relevant part of the master equation for Pt{4>~^, t) is 

tp, = - t^VA^D{ko,t)p, + ... (3) 

lb 
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The volume factor V that appears in the master equation is due to the fact we are 
considering a density matrix which is a functional of two different field configurations, 
(f)^{x) = (j)'^ cos ko.x, which are spread over all space. The time- dependent diffusion co- 
efficient D{ko, t) due to each of the many external environmental fields is then given by 



D^{ko,t) = f dsu{s)[ReGl^{2ko]t-s)+2ReGl^{0;t-s)]. (4) 

J 

In (^), u{s) = cosh^c^o-s when t < tgp, and is an oscillatory function of time when t > t^p. 

Although is oscillatory at all times, for times fit > 1 (until the spinodal one) the 
exponential growth of u{t) enforces a similar behaviour on D^{kQ,t), 

D^{ko,t) ~ (^bJo) exp[2uJot], (5) 

associated with the instability of the k^ mode. For long- wavelength modes D^{t) ~ 
{kBTo/ fiY exp[2fit]. For t > tgp the diffusion coefficient stops growing, and oscillates around 
Dxi^o, t = tsp). 

In our present model the environment fields Xa are not the only decohering agents. The 
environment is also constituted by the short-wavelength modes of the self-interacting field 0. 
Therefore, we split the field as = 0< + 0>, where the system-field 0< contains the unstable 
modes with wavelengths longer than the critical value while the bath or environment- 
field contains the stable modes with wavelengths shorter than fi~^ (in practice, whether the 
separation is made at = /i exactly or at A; ~ /i is immaterial [|15[ by time to, when the 
power of the 0- field fluctuations is peaked at ko fi). This gives an additional contribution 
to the diffusion coefficient. Without the additional powers of A^~^, a high temperature 



resummation of loop diagrams |T6[ is essential to get a reliable 6*++. This is beyond the 
scope of this letter. It will be enough for our purposes to compute an upper bound on the 
decoherence time tr, only considering the external fields Xa- 

We estimate to by considering the approximate solution (|t)f| 



We are following the decoherence time definition of |13| 
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p,[0+ ^p,"[0+ exp\-VT f ds D{ko,s) 

" 

where p" is the solution of the unitary part of the master equation (i.e. without environment) 
and D{kQ^ s) denotes the total diffusion. It is obvious from this (and also from (||)), that the 
diagonal density matrix just evolves like the unitary matrix (the environment has almost no 
effect on the diagonal part of pr). In terms of the dimensionless fields = (0^ + 0;;)/2p, 
and 5 = (0+ - 0-)/2p, we have V = (1/16)^VV^5^- 

The system behaves classically when the non-diagonal elements of the reduced density 
matrix are much smaller than the diagonal ones. We therefore look at the ratio 



Pr[4> + S,4> - 6;t] 



Pr[0,0; t] 



exp 



VT / ds D{ko,s) 



(6) 



p;^[0,0;^] 

It is not possible to obtain an analytic expression for the ratio of unitary density matrices 
that appears in Eq.(|D. The simplest approximation is to neglect the couplings of the system 



field [14]. In this case the unitary density matrix remains Gaussian at all times as 



exp[-^5V'(t)], 
p 



(7) 



p;^[0,0; t] 

where p~^{t), essentially fj? {(j)^)t^ ^ decreases exponentially with time to a value 0{X). A full 
numerical calculation shows that becomes a non-Gaussian function (the associated 
Wigner function becomes non-positive). In any case, in the unitary part of the reduced 
density matrix the non-diagonal terms are not suppressed.^ 

The decoherence time sets the scale after which we have a classical system-field con- 
figuration. According to our previous discussion, it can be defined as the solution to 

1 ^ VT ds D{ko, s) > VT f^" ds DJko, s). 
Jo Jo 



^This should not be confused with the observation that the unitary Gaussian density matrix does 
show classical correlation, whereby the Wigner functional becomes localised in phase space about 
its classical solutions. However, this has nothing to do with eliminating quantum interference 
between different field histories. 
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It corresponds to the time after which we are able to distinguish between two different field 
amplitudes (inside a given volume V). 

Suppose we reduce the couplings ~ A of the system 0-field to its environment. Since, 
as a one-loop construct, F oc (7^ ~ A^, we might expect that, as g, A decrease, then tj) 
increases and the system takes longer to become classical. Although this is the usual result 
for Brownian motion, say 0, it is not simply the case for quantum field theory phase 
transitions. The reason is twofold. Firstly, there is the effect that F oc Tq , and Tq oc A~^ is 
non-perturbatively large for a phase transition. Secondly, because of the non-linear coupling 
to the environment, obligatory for quantum field theory, F oc 0^. The completion of the 
transition finds 0^ ~ r^^ oc A~^ also non-perturbatively large. This suggests that F, and 
hence ^d, can be approximately independent of A. In fact, the situation is a little more 
complicated, but the end result is that does not increase (relative to tgp) as the couplings 
become uniformly weaker. 

In quantifying the decoherence time V is understood as the minimal volume inside which 
there is no possibility for coherent superpositions of macroscopically distinguishable states 
for the field (i.e. there is no 'Schrodinger cat' states inside V). Thus, our choice is that this 
volume factor is 0{fi~^) since fi~^ (the Compton wavelength) sets the smallest scale at which 
we need to look. In particular, characterises the thickness of domain boundaries (walls) 
as the field settles into its ground-state values. Inside this volume, we do not discriminate 
between field amplitudes which differ by 0{^), and therefore we take 6 ~ For we 

set 0^ ~ 0{a/X), where A < a < 1 is to be determined self-consistently from the condition 
that, at time ^d, (0^)* ~ ctv'^- 

Note that the diagonalisation of pt occurs quickly, but not so quickly that /it ^ 1. 
Consequently, in order to evaluate the decoherence time in our model, we have to use 
Eq.(^. We obtain, for the upper bound on t^. 



whereby /xto — ln(?7/TcA/a)). The value of a is determined as a ^ y p/Tc. For comparison, 
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(8) 




we find tsp, for which (0^)* ~ r/^, given by 

2 

exp[2/.t3p] ^ O(^). (9) 

The exponential factor, as always, arises from the growth of the unstable long wavelength 
modes. The factor T^^ comes from the coth.{j3uj /2) factor that encodes the initial Boltzmann 
distribution at temperature Tq ^ Tc. As a result, 

/^isp~ln(-^), (10) 

whereby 1 < fito < fitgp, with 

1 Tc 

/itsp - /it^ ~ -ln(— ) > 1, (11) 
4 /i 

for weak enough coupling, or high enough initial temperatures. Factors of 0{1) have been 



omitted in the argument of the logarithm in (|TT]) and previous equations, requiring that 

This is our main result, that for the physically relevant modes (with small ko) classical 
behaviour has been established before the spinodal time, when the ground states have be- 
came populated. We can say more in that, for an instantaneous quench, nonlinear behaviour 



only becomes important in an interval At, jjAt = 0{1), before the spinodal time [|15[, and 
therefore pr becomes diagonal before non-linear terms could be relevant. In this sense, clas- 
sical behaviour has been achieved before quantum effects could destroy the positivity of the 
Wigner function W^. 
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